Coulomb blockade and transport in a chain of one-dimensional quantum dots 
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A long one-dimensional wire with a finite density of strong random impurities is modelled as a 
chain of weakly coupled quantum dots. At low temperature T and applied voltage V its resistance 
is limited by "breaks": randomly occuring clusters of quantum dots with a special length distri- 
bution pattern that inhibits the transport. Due to the interplay of interaction and disorder effects 
the resistance can exhibit T and V dependences that can be approximated by power laws. The 
corresponding two exponents differ greatly from each other and depend not only on the intrinsic 
electronic parameters but also on the impurity distribution statistics. 



Recently much interest has been attracted by the 
observations of algebraic current-voltage and current- 
temperature dependences in one-dimensional (ID) elec- 



tron systems. The power-laws, I 



V b+1 at high 



V and I oc T a V at low V, appear in a variety of 
conductors including carbon nanotubes 0, 0, H, B, 
nanowires 0, 0, 0, E}- an d polymer nanofibers |lflj . 
This implies that the electron transport in such systems 
is effectively blocked by some barriers at T, V — > 0. It is 
a question of both fundamental and practical interest to 
determine what these barriers are and what controls the 
exponents a and b. 

In the model with a single opaque barrier the power- 
law can indeed appear if electron interactions are strong, 
in which case the effective transparency of the bar- 
rier is governed by the T and V-dependent suppres- 
sion of electron tunneling into a correlated many-body 
state Ull . Similar physics may operate in quasi-lD 
case |l2l Il3l Ihl. Hij . Barring some difference in com- 
putational methods, the existing theories agree that the 
single-barrier exponents should be equal. For an inte- 
rior barrier in a ID wire they are given by a = b = 
2(K~ 1 — 1) where K < 1 is determined by the interac- 
tion strength Experiments with wires a few (im or 
shorter in length Q, 0, U @ support these predictions. 
Yet in wires of length L > 10 /im the ratio a/b is almost 
always larger than unity. It is usually between two and 
five IE S El • To explain this behavior we propose 
that for long wires the single-barrier model may be unre- 
alistic. Indeed, observations of Coulomb blockade (CB) 
at low T [lj, y, llll suggest that a long wire typically con- 
tains multiple impurities that convert it into a chain of 
weakly coupled ID quantum dots. 

We show that for such a system there indeed exists 
a broad parameter regime where the function I{V, T) 
can be approximated by power laws with a 3> b. They 
originate from the statistics of impurity distribution not 




FIG. 1: (a) Illustration of the model. The consecutively num- 
bered dots represent strong random impurities, (b) Typical 
shape of an Ohmic break, (c) A typical non-Ohmic break. 



from Luttinger- liquid effects [ll|. Although in strictly 
ID wires these power laws would be somewhat obscured 
by mesoscopic fluctuations, they may be observable in 
a cleaner manner in quasi- ID systems [HE El El- Un- 
der certain realistic conditions (discussed below) different 
transverse channels of such systems can act as indepen- 
dent ID conductors in parallel, so that our results are 
fully relevant yet mesoscopic fluctuations are averaged 
out. 

The model we study is similar to that of Refs. E3, 
Ihil Hs| and is depicted in Fig. da). Identical impuri- 
ties (e.g., dopants of the same chemical species) are po- 
sitioned along the wire according to a Poissonian distri- 
bution with the mean spacing I 3> 1/n, where n is the 
average electron density. The tunneling transparency of 
each impurity is e~ s <C 1. Electrons occupy a partially 
filled band with a generic (e.g., parabolic) dispersion. 
Their spin is ignored. Electron-electron interactions are 
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short-range, e.g., due to screening by a nearby gate, so 
that the system has a finite (geometry-dependent) ca- 
pacitance per unit length C . The possibility of thermal 
activation and a finite conductance in the CB regime is 
ensured by weak coupling of electrons to the thermal bath 
of phonons. (The conductance depends on this coupling 
through a subleading prefactor, if at all [l^. We ignore 
such prefactors here because they are nonuniversal.) 

The ground state and the low-energy charge excita- 
tions of this model are well approximated by the clas- 
sical limit in which the number of electrons in each 
dot is restricted to be an integer. In the ground state 
it is the integer qj closest to the "background charge" 
Qj = nlj + (£/Aj), where £ is the electrochemical poten- 
tial and lj is the length of the dot. The charging energy 
Aj = (dfi/dn + C _1 )/ij includes contributions both due 
to a finite thermodynamic density of states dn/dfi (the 
analog of the level-spacing in bulk quantum dots) and 
due to electrostatics, as usual. Here /x = £ — U and U 
are the chemical and the electrostatic potentials, respec- 
tively, and the units e = ks = 1 are used. The typical 
charging energy is A = (d/i/dn + C~ 1 )/l. 

Changing the electron number in the dot by m incurs 
the energy cost of E" 1 — (l/2)Aj-m 2 + Ajm(qj — Qj). 
Hence, the equilibrium probability of having qj + m elec- 
trons in j-th dot is given by /™ = exp(— /3i?j™)/Z, where 

f3 = l/T and Z = J2 m ex P(-~P E T^ The dot is in thc 
CB state if \q 3 - Qj\ - 1/2 > T/Aj. In this case only 
m = 0, ±1 are important. As common in CB literature, 
one can view = fj, ± E^ 1 as the energy of the lowest 
unoccupied (highest occupied) "single-particle" orbital. 
Important in the following is the probability P(e + ,e~) 
of having no charge excitations within a given energy 
interval, E^ 1 > [2(|. Function P depends only on 
e = e + + e~ and can be easily calculated: 

PM !_«(!_.-*/.) (a) 

(the integration is constrained to Aj > e). It behaves as 
A/(2e) at e > A, and so the average charge excitation 
gap (e) = —JedP(e) is logarithmically divergent |2l|. Be- 
low we turn to the study of transport and show that this 
divergence and the constrains of the ID geometry lead to 
new conductivity laws. 

We focus on the case of low voltages, V < V c = 
LA/ (21), where the wire is in the CB state j2^. Electron 
transitions between any two dots j and k occur by the 
combination of quantum tunneling and thermal activa- 
tion. For \j—k\ = 1 the tunneling is said to be sequential; 
otherwise, it is referred to as cotunnelin g |1 4| . In prior 
studies of random quantum dot arrays [13, 12a, |2g, L=JZJ| 
the analogy between this transport mechanism and the 
variable-range hopping |U 0] (VRH) has been noted. 
Thus, the typical length of a hop \j — k\ is established 
from the competition between activation and tunneling. 
Longer hops allow electrons to find transitions of lower 



activation energy but involve a higher tunneling action 
s\j — k\. The crucial point is then as follows. It is known 
that ID VRH (of non-interacting electrons) is special: 
the conductance is limited not by typical hops but by 
rare highly resistive spots — "breaks" [2^, |2?J HO, EhJ- 
Therefore, one may anticipate that some type of breaks 
operate in our model as well. Indeed, the calculation be- 
low shows that at not too low T the role of breaks is 
played by rare clusters of densely spaced impurities, see 
Fig. Ha) 0. 

The calculation closely parallels those of Refs. [30( and 
[3l| . It starts by noting that if breaks exist, they act 
as transport bottlenecks where the most of the applied 
voltage drops. In contrast, the dots between two adja- 
cent breaks are in a quasi-equilibrium. Their electro- 
chemical potentials £j are nearly equal while the occu- 
pation factors of their charge excitations are given by 
/• 1 ~ exp[±/3(/Uj — where fj,j is the local chemi- 

cal potential. This allows one to get a simple formula 
for the net current Ijk from a dot j to another dot k. 
For the most relevant case of a large activation energy, 
Ejk S> T, it coincides with the standard expression of 
the VRH theory: 

I jh ~ I e- S ^e-^ sinh/3(0 - (2) 

E kj = \ min (|ej + \e\ - fi k \ + |ej - e£|), (3) 

where Iq = const is some prefactor 33] . Consider first 
thc Ohmic regime. Here fij ~ fi = const and the effective 
resistance Rj k = (£j —£,k)/Ijk f° r each link of the hopping 
network can be defined: Rj k = Roexp(f3Ejk + s\j — k\) 
with i? = T/I Q . Denote by P u (u) the probability per 
unit length of finding a break of resistance at least Roe u , 
i.e., such a configuration of dots that Rjk > Roe u regard- 
less of how j to the left of the break and k to the right of it 
are chosen. For u s this break contains many, at least 
u/s, dots. We can then talk about its shape, i.e., smooth 
envelope functions e ± (z) > Osuchthat ±(e^— fi) > e ± (j) 
for all j in the break, see Fig. ^b) . The probability in 
question can be approximated as follows: 

P u (u) ~ maxJ]P[ £ +(i) +E~(j)], (4) 

3 

where the maximum is sought under the constraint Rj k > 
Roe u for all (i, k) pairs. Verbatim repetition of the steps 
taken in Ref. |30| establishes that the break has a familiar 
diamond-like shape H3,|^|, see Fig.^b). 

It requires some effort to finish the derivation of P u for 
an arbitrary u. However, for Tu ^> A the final result can 
be obtained if the actual shape of the break is approxi- 
mated by a rectangle of the same width and height, i.e., 
if we consider a break of ./V = u/s dots with the charging 
energy Tu each. This gives P u (u) ~ [P(Tu)] N , and so 

LnP„(u) ~ -(u/s)bx(Tu/A), Tu > A. (5) 
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A non-Ohmic break, where Vb = £j — £fc ^ T for 
j and k on the opposite sides of the break, can be 
treated similarly. The necessary steps are spelled out in 
Ref. [3l| • Namely, consider a state of a given fixed cur- 
rent I = ioexp(— uj), so that the free parameter is the 
voltage Vb generated by a break. Denote by PviVb) the 
probability per unit length of finding a non-Ohmic break 
with the voltage drop of at least Vb- Once again approxi- 
mating the break shape by a rectangle e ± = Tuj + (V>/ C) 
with C ~ 1 and using Eq. we find this time 



\nP v (V b ) ~ \nP u { Ul ) - Vb/{sCT). 



(6) 



This entails that the average voltage drop V» across a 
non-Ohmic break is V* = s£T. Thus, if V S> V*, many 
breaks must contribute. The derivation of I(V) is espe- 
cially simple in this case, so we prefer to finish it first 
and turn to the Ohmic one later. 

The total voltage V is the sum of voltages generated 
by all non-Ohmic breaks (contribution of Ohmic ones is 
subleading). For V S> V* (many breaks) the sum can 
be replaced by the integral V = L J dVbPv(Vb). Using 
Eqs. J5J) and Jfjjl we get a transcendental equation for ui 

ui = sln( Cl LT/Vl)/ln(T Ul /A), 7* < V<€ V c , (7) 

where ci is a slow dimensionless function of V and T. We 
conclude that the dependence of I on V is close to the 
power law with a weakly T and 7-dependent exponent 



b ~ s/ln(Tu//A), T>A/s. 



(8) 



This is the first of our two main results. As V decreases 
to a number ca. V, the resistance of a wire becomes 
dominated by the single largest break, and so it exhibits 
strong ensemble fluctuations. To illustrate how to handle 
this case, let us consider next the Ohmic regime where 
such fluctuations are the strongest. 

In this regime the resistance R of the wire is the sum 
of Ohmic resistances Rj — i?o exp(itj), where Uj's form 
a set of Nb ~ L/lb independent random numbers with 
the probability distribution function (PDF) lbP u (u) each 
[Eq. ©]. Parameter l^ 1 = fP u (u)du enters the final re- 
sult only under the logarithm, so it need not be specified 
precisely. 

Following j^II, our goal is to compute Pr, the PDF of 
R. We begin by noting that the PDF p of Rj behaves as 



P (Rj) ~ r: 



-l-a 



a = s- 1 ln[(T/A)ln(R j /R )], (9) 



cf. Eq. J5j. Since a(Rj) is a very slow function, p(Rj) 
is basically a power-law. Hence, the random variable 
R > is a Levy random walk. A celebrated theorem 
of the probability theory immediately tells us that in 
the limit of large Nb the PDF of R approaches the sta- 
ble Levy distribution pjjf (often encountered in plasma, 



astro, and atomic physics, as well as in biology, eco- 
nomics, and reliability theory). For a <C 1 it becomes 
identical to Frechet extreme- value statistics Pr(R) oc 
fl-i-a, exp[-(R*/R) a *]. In turn, the PDF of the con- 
ductance G — 1/R has the form 



P G {G) = const x G Q *- 1 exp[-(i?,G) Q *], 



(10) 



so that its first two moments are given by 

{G)R* =T(1 + a- 1 ), (G 2 }Rl=T(l + 2a: 1 ), (11) 

T(z) being the Euler gamma-function. 

The equations for a* and i?* are obtained by the stan- 
dard procedure [35j (see also Ref. |30j'): 

a, ln(i?„/i? ) = ln(AT 6 /a»), a, = a(R*). (12) 
Combined with Eq. 0, they yield 



\n(RJR ) 



s ln(L/a*4) 



ln[(Ts/A) ln(L/aJ b )] ' 



(13) 



Not surprisinglvjthe form of Eq. (jlOJI is identical to that 
found in Ref. 30] for another model; however, the In- 
dependence of i?* in that case is different from Eq. 113fl . 
If only a limited range of T is accessible, which is often 
the case in the experiment, the latter can be approxi- 
mated by a power law with a weakly varying exponent. 
In turn, this entails (G) oc T a , where a is given by 



d\n(G) ^ \n(c 2 L/l) 
d\nT ~ ln(Ts/A) 



(14) 



and C2 is an algebraic function of T [361 ]. This is our 
second main result. As emphasized above, a»i)> 1. 

Equation (JT3J) applies at T > A/[s ln(L/Z)] and for 
wires of length I <C L <C Zexp(e s ). The lower bound 
guarantees that the wire contains a large number of im- 
purities. The upper bound, which is likely to be satisfied 
in practice, ensures a* -C 1, so that Eq. l(TU|) holds. 

The PDF of Eq. ltTU|) describes the statistics of G in an 
ensemble of wires at a given T. However, fluctuations of 
comparable ma gnit ude would appear in the same wire as 
T varies They would be superimposed on the 

backbone dependence G a T°. According to Eq. (fTTfl . 
the amplitude of such fluctuations is large in the limit 
we considered, a* <C 1. Let us however estimate how 
small a* can really be. For the experimentally accessible 
regime of \n(R*/R) < 10 and Nb, L/l of a few hundred, 
Eq. I|12|) gives a* ~ 0.5. In this case, the fluctuations of G 
are still significant [cf. Eq. Qllfl]. and so verification of our 
predictions in truly ID systems would require statistical 
analysis of the data j3]} and/or ensemble averaging. 

At this point we draw attention to the fact that the 
cited experiments 0, H, 0, El were actually performed 
on quasi-lD systems with a large number N of trans- 
verse channels. Suppose that conductances of different 



4 



channels are additive and statistically independent, then 
for a* = 0.5 the relative fluctuations of the total G will 
be y/b/N. This amounts to 10% for a representative 
number of iV = 500, and so the power laws should no 
longer be obscured by fluctuations. The above scenario 
certainly needs an experimental verification. It will not 
hold for materials where transverse channels are coupled. 
However, there are systems in which it is reasonable, e.g., 
those with a large concentration of neutral (in general, 
short-range) impurities and strong anisotropy of longitu- 
dinal and transverse transport, see Ref. 17]. 

Finally, let us comment on the expected values of a, 
b, and s. For a ID wire the result s = — Int + {K^ 1 — 
1) ]n(E%/5jSj-i) can be obtained HIEU^. Here t < 1 
is the bare tunneling transparency of the impurity and 
Ef, Sj are, respectively, the high and the low energy 
cutoffs for the Luttinger-liquid effects. The latter cutoff 
Sj = KAj oc is the energy spacing of the collective 
neutral excitations of the jth dot. In comparison, for an 
impurity in an infinite wire, the low-energy cutoff is set 
by either V or T 0] . Since s depends on lj 's only loga- 
rithmically, our approximation that s is the same for all 
impurities is justified. Note that Sj ~ KTu oc T\n{L/l) 
for important Ohmic breaks, so that, nominally, s has the 
same logarithmic T-dependence as in the single-barrier 
case (llj. Nevertheless, due to the (modestly) large log- 
arithmic factors h\(Ep/Sj) and h\(L/l) both predicted 
exponents a and b would noticeably exceed the single- 
impurity one, in this case, 2(K~ 1 — 1). 

Semiconductor nanowires, polymer nanofibers, and 
carbon nanotubes appear promising for a continuing in- 
vestigation of the physics we discussed. We call for a 
systematic study of transport in such wires. 
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